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Abstract 

Some pressure and pressure-velocity correlations in a direct numeri- 
cal simulations of a three-dimensional turbulent flow at moderate Reynolds 
numbers have been analyzed. We have identified a set of pressure- 
velocity correlations which posseses a good scaling behaviour. Such a 
class of pressure-velocity correlations are determined by looking at the 
energy-balance across any sub-volume of the fiow. 

According to our analysis, pressure scaling is determined by the dimen- 
sional assumption that pressure behaves as a "velocity squared" , unless 
finite-Reynolds effects are overwhelming. The S0(3) decompositions of 
pressure structure functions has also been applied in order to investigate 
anisotropic effects on the pressure scaling. 

1 Introduction 

Scaling in turbulent flows is one of the most challenging open issue in fluid- 
dynamics |T|]. Typical problems concern both the understanding of the ideal 
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case of isotropic and homogeneous turbulence in the hmit of high Reynolds 
numbers [0, 0] or more realistic and applied situations with anisotropic and 
inhomogeneous statistics (for recent examples see In 1941, Kolmogorov, 

used a clever application of dimensional analysis to predict that the scaling of 
velocity increments in the inertial range should have a power law behaviour 
depending only on the averaged energy dissipation in the flow, e. Namely, for 
structure function of order p, i.e. the p-th moment of a velocity difference 
across a distance R we have: 

Sg{R) = {{v{x + R)- v{x)y) ~ 6"/=^ R''^^ (1) 

with 7] <^ r Lq where rj is the dissipative scale and Lq is the typical external 
scale where forcing acts. Let us notice that in (|ID we have explicitly neglected 
any tensorial structure in the velocity field such as to stress the typical dimen- 
sional character of the Kolmogorov theory, i.e. the scaling properties must be 
the same for any observable which has the same physical dimension and which 
is built in terms of local field increments, 5rv{x) = f (x + r) — v{x). 
Kolmogorov theory, as previously summarized, is quantitatively wrong. Ex- 
periments and numerical simulations show a quantitative disagreement with 
the dimensional prediction p/3 for the scaling exponents. For example, the 
longitudinal velocity structure functions: 

SI (r) = ( {v, (x + R) - V, (x)) 4 ') ~ (2) 



show a power law behaviour with a set of exponents C,v{p) non linear in p. 
The failure of the dimensional estimate p/3 goes under the name of anomalous 
scaling. 

Many problems naturally arises as a consequence of the failure of the main 
Kolmogorov prediction. The main open problem is to find an analytical way 
to calculate from first principle the anomalous exponents, a problem which is 
still out of control but for the case of anomalous exponents characterizing the 
statistics of passive quantities advected by gaussian velocity fields ff\-^ ■ 



Another interesting question, opened by the failure of Kolmogorov dimensional 
prediction, consists in the possibility that local observable with the same phys- 
ical dimensions but with different tensorial structures have different scaling 
properties. About this point, there are some experimental and numerical ev- 
idences regarding different possible anomalous behaviour of longitudinal and 



transversal structure functions |]13[ even in isotropic turbulence. Somehow re- 
lated to this issue is also the apparent different anomalous scaling between 
the coarse grained averages of dissipative quantities like enstrophy and energy 
dissipation WM. On the other hand, on the basis of a S0(3) decomposition of 
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velocity correlation functions, other authors [jT5[ claim that the supposed dif- 
ferent scaling of quantities like transversal and longitudinal structure functions 
can only be due to spurious sub-leading non-isotropic effects, i.e. in isotropic 
high- Reynolds numbers all components of the same tensorial observable should 
have the same -maybe anomalous- scaling behaviour. A first numerical support 
to this claim has been presented in the analysis of a channel flow simulation in 

i- 

Even more complex is the situation when multi-point pressure correlations is 
involved P,[l3,|T^. Dimensionally speaking pressure is just a velocity squared, 
and Kolmogorov-like argument can be easily generalized to the case of pressure 
structure functions, Fq{r). Indeed, a simple application of dimensional analysis 
leads to [0: 

F,(r) = (|P(x) - P(x + r)n ~ e^^? (3) 

where as usual, all distances are supposed to belong to the inertial range of 
scales. Of course, intermittency will also affect pressure scaling. By follow- 
ing the straightforward hypothesis that pressure can be treated as a velocity 
squared one would be tempted to assign the same intermittency exponents of 
the velocity field to the pressure scaling, i.e. to replace with 

This prediction is just a simple consequences of the assumptions that all ve- 
locity correlations have the same scaling behaviour supposed that all distances 
involved are in the inertial range and that the statistics is locally isotropic. 
Such a dimensional ansatz has been questioned on the basis of a phenomeno- 
logical argument in [^, some numerical support to the latter argument have 
been recently presented in ||18|| . 

In this paper we will mainly present some numerical evidences that indeed the 
dimensional ansatz (HI) is correct. It is well known that this must be the case 



at least for for g = 2 in (^). In this case, there exist an exact relation pO | 
which connect the scaling of the second order pressure structure function with 
a linear integral combination of fourth-order velocity structure functions. The 
problem is if the exact result can be simply extrapolated to other pressure- 
dependent observable and, in the case, how strong finite-Reynolds effects can 
be. Indeed, one may argue that pressure feels strongly non-local effects, being 
just the inversion of the Poisson problem AP = —didjViVj, and therefore the 
assumptions of independence from large scales and/or from boundary condi- 
tions may not be satisfied even at very high Reynolds numbers. Indeed, to 
our knowledge, neither experimental studies nor numerical simulations have 
ever been able to make a firm quantitative statement about pressure scaling 
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properties |16|-[T8|, . In this paper we show that it is possible to find a set of 
velocity-pressure observable which have indeed a quite good scaling behaviour 
in agreement with the dimensional ansatz (^) also at moderate Reynolds num- 
bers. 

Scaling in turbulence is particularly difficult to test in both experiments and 
numerical simulations. Experiments reaches high Reynolds numbers by paying 
the price to have a very limited set of information on the whole velocity fields, 
typically only a long time series of one velocity components in a few spatial 
points. Moreover, in most cases, there is not a precise control of the degree 
of isotropy and homogeneity in the fiow. On the other hand, numerical simu- 
lations have a perfectly controlled lay out, the velocity field is exactly known 
at any point, but the maximum reachable Reynolds number is still order of 
magnitude smaller then in typical experiments ||25[ . 



Nevertheless, numerical simulations, if exploited in a clever way, are the only 
tool where complex measurements can be performed. Therefore, questions like 
the dependency of scaling properties from the tensorial nature of the observ- 
able can, up to know, be investigated only in numerical data base. 
In this paper, we present a detailed analysis of pressure scaling and pressure- 
velocity correlations scaling in a set of moderate Reynolds number simulations. 
Starting from the analysis of the energy transfer in real space we propose a set of 
pressure- velocity observable which show better scaling properties then the usual 
pressure structure functions. We presents quantitative evidences that indeed, 
while pressure structure functions are strongly affected from Reynolds numbers 
effects, the pressure- velocity correlation functions we investigated have a fairly 
good scaling behaviour, even at modest Reynolds numbers, in agreement with 
the hypothesis that pressure "behaves" like a velocity squared. In order to 
understand whether the bad scaling behaviour detected in the pure-pressure 
structure functions is due to spurious anisotropic sub-leading effects we also 
present some results on the SO (3) decomposition of the pressure field. 
The paper is organized as follows. In section 2 we summarize the known anal- 
itycal result which connect the second order pressure structure function to the 
integral linear combination of fourth-order velocity correlations and the exper- 
imental and numerical attempts to test the relation. In section 3 we introduce 
the set of pressure-velocity correlations which should have a better scaling 
properties on the basis of a simple argument based on the energy transfer of 
Navier-Stokes equations in the real space. In section 4 we present the analysis 
of our numerical data base. In section 5 we briefiy comment on the analysis of 
non isotropic fiuctuations. Conclusions follow in section 6. 
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2 Pressure structure functions 



Under the assumptions of local isotropy, local homogeneity, incompressibility, 
and by use of Navier-Stokes equation, one can relate the second order pressure 
structure functions, F2{r), to some fourth-order velocity structure functions 
ig. Namely: 



1 4 r'^ 

F2(r) = --Dniiir) + -r' y-^ [Duii{y) + Dp^p^iy) - 6A177] % 



4 /"CO 

+ 2 J [Dpi3pp{y) - 3Dii^^] dy (5) 

where the fourth-order structure function is 

D^jklif) = {{ui- u[){uj - u'.){uk - -u'l)) 

and where for simplicity we have used primed variables to express velocities at 
the position x' = x + r and where i,j, k,l isl if the velocity component is paral- 
lel to the separation vector, r, and 2, 3 otherwise. Subscripts /?, 7 denote either 
2 or 3. Of course, implies that whenever the fourth-order structure func- 
tions entering in the above expressions are all dominated by the inertial-range 
intermittent scaling behaviour, Dj j /(. ;(r) ~ r'^^^^\ then also the second order 
pressure structure functions should scale with the exponent Cp(2) = Ci;(4)- Re- 
lation (j^) have been carefully tested in numerical simulations without any ap- 
preciable deviations [16, 1^. Nevertheless, the overall scaling behaviour of the 



pressure structure function is very poor. Similarly, the analysis of experimen- 
tal data fl^ does not show any power law behaviour for the pressure structure 
functions even if the Reynolds number was extremely high [Rex ~ 10000). In 
the latter case, authors tried to explain the difference between pressure scaling 
quality and velocity scaling quality by invoking a possible different scaling for 
the different velocity correlations entering in the RHS of (|^), leading to the 
final prediction that pressure structure functions is made in terms of different 
power law contributions with slightly different exponents. The resulting super- 
position of power laws would be the responsible of the poor observed scaling 
behaviour. This statement would anyhow contradict the theoretical prediction 
made in terms of the SO (3) decomposition which forbids different component 
of the same tensorial observable to scale differently in a isotropic ensemble. 
Another interesting remark consists in the strong cancellation among the dif- 



ferent contribution of (^) observed in numerical simulations ||16|,p!7|: the LHS 
of is more then an order of magnitude smaller than the single different 
contributions entering in the RHS. One, cannot exclude apriori the possibility 
that there exist an almost perfect cancellations of all leading scaling terms of 
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all contribution appearing in the RHS of (^), even if such a perfect cancellation 
would call for some unknown physical interpretation. 

More probably, the cancellation is not perfect but strong enough to hide 
completely the pressure scaling at the available experimental and numerical 
Reynolds numbers. 

On the other hand, the possibility that pressure-increments behave as velocity- 
increments, 5P ~ 5f , instead than as a velocity-increment squared has been 
recently proposed |]T8[. This would violate the exact results previously reported 



and therefore cannot be correct unless strong anisotropic effects are present at 
all scales. 

In order to better assess the pressure statistical properties we present in the 
following, some results for pressure-dependent observable. This observable do 
not posses the strong cancellation properties showed by structure function. 



3 Pressure-velocity correlation 

Let us start by looking at the energy balance inside any volume V of the flow. 
From the Navier- Stokes eqs we obviously have: 

dtEy + J fi(x)t>j(x)9jf j(x)(ix + J Vi{x)diP{x)dx = u J t>j(x) At>j(x)(ix (6) 

where with Ey = 1/2 Jy dx.ViVi we denote the total energy in the sub- volume 
V. Let us notice that the two terms in the LHS of (H) can be obviously written 
as the fluxes across the boundaries of V by using Gauss theorem: 

/ fi(x)t;j(x)9jfj(x)rfx = / dJ:njVj{x.)v'^{x) = ^Q^{vv^) (7) 
f v,{^)diP{^)d^^ f rfSn,t;,(x)P(x) = <l>n^(vP) (8) 

JV Jftv 

where with Ui we denote the versor perpendicular to the infinitesimal surface 
on the boundaries of V. Rewritten in this way, relation is just a simple 
restatement of the conservation of energy: the total energy change inside a 
volume is given by the flux across the volume surface and by the energy dissi- 
pation inside the volume. Let us now use this simple fact in order to extract 
some useful guess about scaling properties of velocity-pressure correlations. 
Let us consider a very particular class of volume V, i.e. a cylinder with an 
infinitesimal squared basis of surface e and with a finite axis in the direction 
of R. In the limit when the basis becomes smaller and smaller, the flux across 
the lateral sides goes to zero because contributions from two opposite walls are 
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equal but with different sign. The only contributions to the total ffux come 
from the two infinitesimal basis and can be written as 

^ny{y^v^) = e{vi{^y{^) - t;,(x + R)t;2(x + R))4 (9) 

for the flux involving the velocity correlation and as 

$f,,.(vP) = e(t;,(x)P(x)-t;,(x + R)P(x + R))i?i (10) 

for the flux involving pressure-velocity correlations. In both cases we have 
exploited the fact that the two infinitesimal basis are centered in x and in 
X + R, i.e. their versor is oriented along R. Similarly the two volume integral 
giving the time variation of the total energy and the energy dissipation becomes 
two linear integral times the infinitesimal basis area e: dt{e J-^dsv'^{x.)) and 
(e JYidsVi{x.)Avi{x.)) respectively, where with ds we parametrized the segment 
going from x to x + R. 

Let us now assume that all the four observable entering in the energy balance 
have the same statistical behaviour. This is somehow a "local Kolmogorov 
refined hypothesis": we link the scaling of the local energy dissipation to the 
scaling of some particular third order velocity correlation and to the scaling 
of a velocity-pressure correlation. The claim is therefore that the particular 
structure functions emerging from our flux analysis should scale exactly like the 
coarse grained energy dissipation, i.e. should have an anomalous scaling like the 
usual longitudinal structure functions which satisfy the original Kolmogorov 
Refined Hypothesis. In the next section we present some numerical data in 
support of this claim. 



4 Numerical Analysis 



The data set we are going to analyze has been obtained from a direct numerical 
integration of NS eqs using a pseudo-spectral method with dealiasing on a grid 
of 128^ points. The forcing was implemented isotropically on all wavevectors 
with |fc| < 1 such as to enforce the k~^^^ spectrum at small wavectors We 
have analyzed about 100 configurations stored each eddy-turn over time. The 
simulation has a Re\ = 70. 
Let us denote with: 



S7 (R) 



Sf{R) 



Vi{x)v'^{x) -Vi{x + R)t;2(x + R)) P, 
t;i(x)P(x) -Vi{^ + R)P(x + R)) R 



q/3 
9/3 



:iii 
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the two different structure functions wliich can be made in terms of tlie two 
flux quantities defined in tlie previous section. Let us notice tliat in ([TT|) botli 
^™ (i?) and Sf{R) liave been defined as tlie q/3 power of tlie original fluxes 
such as to have the same dimensions of S'^{R). 

Let us start by showing in Fig.l the strong cancellation effects present in the 
pressure structure functions Fg{R) with respect to the velocity longitudinal 
structure functions with the same physical dimensions S2q{R)- In Fig. 1 we 
show the log-log plot of Fi{R) and of 5*2 (-R), as one can see the overall ampli- 
tude of pressure fluctuations is about one order of magnitude smaller than the 
velocity fluctuations. This is just to confirm that pressure by itself is a much 
weaker signal than the usual velocity correlations. 

As one can see in Fig. 1 the scaling is quite poor, as one can expect in any 
DNS. As usual, in order to extract quantitative statement about scaling expo- 
nents one has to exploit the Extended Self Similarity (ESS) property enjoyed 
by homogeneous and isotropic turbulent flows [^. ESS consists in looking for 
relative scaling of two different observable. Usually, one takes two structure 
functions of two different orders, i.e. in the case of longitudinal structure func- 

tions s:;{R) \s:;,{R)f'^'^^'^''^''\ 

Let us now define the same relative scaling for the two generalized structure 
functions defined in (pT|). 



S^/{R) ~ [S^/{R)\ , ~ [Sg^{R)\ ^"^'''^ (12) 

In Figs. 2,3,4 we show the ESS plot for the two generalized structure functions 
and for the usual longitudinal structure functions respectively, with q = 1, q' = 
2. As one can see the all ESS plots showed a scaling behaviour consistent with 
the usual homogeneous and isotropic high Reynolds value which give for the 
relative exponents: Ci)(2)/Cd(1) = 1.92±0.02 [^. Similar agreements are found 
for higher order moments (not showed). 

The scaling ansatz assumed for the generalized structure functions seems there- 
fore quite well satisfied. These findings support the fact that pressure does not 
behave abnormally as far as its "dimensional" scaling properties are concerned. 
Indeed, pressure-velocity correlations behaves exactly like velocity-velocity cor- 
relations once pressure is counted as a "velocity squared". Nevertheless, The 
only realistic way to perform a quantitative statement about scaling exponents 
is to study the logarithmic local slopes of (|1^). Only when logarithmic local 
slopes show a fairly constant behaviour one can really speak about scaling. In 
Fig. 5 we show the logarithmic local slopes of (|T2|) for g = 1, g' = 2 together 
with the corresponding quantities measure for the longitudinal velocity struc- 
ture functions. In order to show that the hypothesis that pressure-increments 



„2 , 
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behave as a linear velocity- increment, 6P ~ 6v, is definitely ruled out by our 
data we also show in Fig. 5 the logarithmic local slope of the ESS applied to 
pressure structure-functions for q = 1, q' = 2. As one can see, while the three 
slopes measured on the flux structure functions and on the longitudinal struc- 
ture functions agree perfectly with the high-Reynolds numbers measurements, 
the pure-pressure structure functions is definitely much poorer. In Fig. 6 we 
plot the same of Fig. 5 but for a different choice of moments, q = 2, g' = 4, 
for both fluxes and longitudinal structure functions and with g = 1, g' = 2 for 
the pure pressure structure functions. In this way we are comparing quanti- 
ties with exactly the same dimensional properties. Again, while the fiux-made 
structure functions, Sg^{R), S"™ (i?) and the longitudinal structure functions, 
S^{R) have the same local slope the pure-pressure result obtained on the ESS 
of Fq{R) shows a poorer and different scaling. 

A few comments are now in order. On one hand, we see from Figs. 5 and 
6 that the simple local-refined Kolmogorov hypothesis derived in the previ- 
ous section is correct, i.e. fluxes ([TTl) have the same scaling properties of the 
usual longitudinal structure function in homogeneous and isotropic turbulence, 
confirming that these observable with the same physical dimensions and built 
in terms of local field increments scale in the same way. On the other hand, 
(see Fig. 6) pure-pressure structure functions seem to violate the previous 
statement despite of the fact that in this case there even exist an exact result 
(^) supporting it. Why pure-pressure correlations show this strong deviation 
from the straightforward dimensional estimate? One possible explanation is 
connected to the -possible- lack of isotropy in the statistics. Any isotropically 
forced DNS is affected by possible non-isotropic fluctuations induced by the 
discretization of the numerical grid. In the following section we have analyzed 
non-isotropic effects on both velocity and pressure fluctuations. 



5 Anisotropic effects 

The exact relation which connect the second order pressure structure function 
with a linear integral combination of fourth order velocity structure function 
(^) is correct only in the isotropic and homogeneous case. In order to test 
the degree of isotropy of our simulation we have proceeded in a systematic 
decomposition in terms of the irreducible representations of the SO (3) symme- 
try group [|, |, The S0(3) decomposition is particularly simple to apply 
to scalar observable, i.e. observable whit all vectorial indexes contracted, like 
pressure structure functions or longitudinal structure functions on the kind an- 
alyzed in this work. In these case, the SO (3) decomposition is nothing but a 
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decomposition in spherical harmonics. For example, the longitudinal structure 
functions, 3^(11) =< (fi(x + R) — t>j(x))_Rj)'^ >, can be decomposed as: 

s:;{K)^Y.srimU^) (13) 

jm 

where now, we have explicitly considered the possibility that the undecom- 
posed structure functions depend on the whole vector R and not only on its 
magnitude as in the previous sections when isotropy was assumed. The coeffi- 
cient of the decomposition, S^"^{R) depend only on the magnitude of R and on 
the two "quantum" numbers j, m which labels the properties under rotations of 
the Yjm eigenf unction. Obviously, in the case of perfect isotropy we would have 
only one projection alive, i.e. the projection on Yqo- The SO (3) decomposition 
here summarized as been already used in some experimental and numerical 
data analysis to properly disentangle the anisotropic effects from the isotropic 
ones [^, ^, |T^, |T3[ . In our case, the relative amplitudes of S'™(|i?|) with respect 
to the anisotropic fluctuations S^g"^{\R\) with j > 0, —j < m < j , gives a direct 
quantitative estimate of the degree of anisotropic fluctuations for any scale \R\. 
In Fig. 7 we show the log-log plot of the undecomposed second order longitudi- 
nal structure functions and of its projection on the fully-isotropic eigenfunction 
Foo- As it is possible to see, despite of the isotropic forcing used in the sim- 
ulation, the flnite-size effects introduced by our computational grid are quite 
important at large scales: the projection, 5'2°(|-R|), shows a deflnitely better 
scaling than the undecomposed structure function already as a functions of 
the real separation R, i.e. without using ESS. This dramatic effect was already 
observed in a similar application to the decomposition of velocity fluctuations 
inside a channel |Q. Fig. 7 deflnitely show that the SO (3) decomposition can 
help in cleaning scaling properties also in "quasi-isotropic" simulations. 
On the other hand, the situation is quite different when the same decomposition 
is applied to the pressure structure functions: 

F:;{K)^Y.Frimu^) (14) 

jm 

Let us notice that pressure is a quasi-isotropic observable also for strong anisotropic 
velocity conflguration. Indeed, being the solutions of a Poisson problem, pres- 
sure is always an average of velocity fluctuations on all spatial directions. This 
simple considerations is perfectly verifled on our numerical simulation. In Fig. 
8 we show the undecomposed second-order pressure structure function together 
with its projection on the fully isotropic harmonics, One can hardly 

detect any differences, suggesting tha anisotropic fluctuations cannot be the 
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responsible for the poor scaling observed in the previous sections for the pure- 
pressure structure functions. 

The only possibility to reconcile the exact result with the poor scaling 
agreement between pressure structure functions and velocity structure func- 
tions is, in our opinion, to invoke strong Reynolds effects. Indeed, the reso- 
lution of the Poisson problem certainly introduces strong non-local effects on 
the statistics. Non-locality may also translate in strong long-range effects in 
the Fourier space as far as the importance of boundary conditions and forcing 
on the inertial range properties are concerned. If this is correct, there are not 
reason to expect good scaling properties for pure-pressure correlations, unless 
Reynolds number is high enough to recover also in laboratory experiments an 
almost- "infinite" inertial range extension. 



6 Conclusions 

We have analyzed some pressure and pressure-velocity correlations in a Direct 
Numerical Simulations at moderate Reynolds numbers. We have derived on the 
basis of a simple analysis of energy transfer properties across any sub-volume 
in the real space what we call a "local" -Refined Kolmogorov Hypothesis. We 
have identified a set of pressure- velocity correlations which should have a good 
scaling behaviour because connected via the local-RKH to the scaling of the 
energy dissipation coarse grained on inertial range scales. 
We have showed that our scaling hypothesis is well verified, while pure pressure 
correlations feels strong Reynolds effects. According to our analysis pressure 
scaling is perfectly determined by the dimensional assumption that pressure 
behaves as a "velocity squared", unless the finite- Reynolds effects are over- 
whelming. We do not find any sign which could support the fact that pressure 



differences behave as velocity differences as proposed in . 
We have also applied the SO (3) decompositions to the pressure structure func- 
tions in order to show that poor scaling properties showed by pure-pressure 
structure functions are not connected to anisotropic fluctuations. 
We acknowledge some help by I. Mazzitelli in the SO (3) analysis. LB and 
FT have been partially supported by INFM (PRA-TURBO) and by the EU 
contract FMRX CT98-0175, 
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Figure 1: Log-log plot of Fi{R) +, and of S2{R), x. Notice that the pressure 
structure function is about one order of magnitude smaller than the velocity 
structure function at large scale. 
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Figure 2: ESS log-log plot of 82'"^ (R) vs Sl^^{R), superimposed is the straight 
line with the isotropic and homogeneous high- Reynolds slope Q/Ci — 1-92. 
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Figure 3: The same as in Fig. 2 but for S'^{K) vs S'^{K). 
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Figure 4: The same as in Fig. 2 but for the longitudinal structure function, 
Sl{R) vs Sl{R). 
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Figure 5: Logarithmic local slope of: (*) S^^^R) vs SJ^R)- (□) S^^{R) vs 



SfiR); (x) S^^{R) vs for g = 2 and q' = 1. 



Notice that all the above 



values are in perfect agreement with the high- Reynolds number value 1.92 
(straight line), while the logarithmic local slope for the pure-pressure structure 
functions F2{R) vs -Fi(-R), (+) is different. The error bars are estimated by 
looking at the fluctuations over the first half and the second half of the whole 
set of configurations. 
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Figure 6: Logarithmic local slopes of: (□) ^^'(i?) vs S^fiR); (*) Sf{R) 
vs SfiR); (x) S"g{R) vs for g = 4 and q' = 2. Notice that, as in 

Fig. 5, the flux-based structure functions have the same scaling behaviour of 
the longitudinal structure function in agreement with the high- Reynolds regime 
(straight line). Here we present a comparison with the Pressure structure func- 
tion with the same physical dimensions of the flux-based structure functions, 
i.e. F2{R) vs Fi{r), {+), still the pure-pressure structure function seems to 
have a different local slope. 
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Figure 7: Undecomposed second order velocity structure functions >S'2(R) mea- 
sured on the plane x — y, (x); and the projection, (R) (+), on the isotropic 
eigenfunction. The straight line has the high-Reynolds slope Q — Notice 
that already in the i?-space, the SO (3) decomposition improve the overall scal- 
ing behaviour. The two curves have been shifted along the y axis for the sake 
of presentation. 
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Figure 8: Undecomposed second order pressure structure functions F2(R) mea- 
sured on the plane x — y, (x); and the projection, F2^{R) (+), on the isotropic 
eigenf unction. Notice that, at difference from Fig. 7, here the decomposed 
and undecomposed structure functions are almost identical, indicating a high 
degree of isotropization in the pressure statistics. 
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